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Abstract 

We discuss the generation of the baryon asymmetry by a strong first order elec- 
troweak phase transition in the early universe, particularly in the context of the 
MSSM. This requires a thorough numerical treatment of the bubble wall profile 
in the case of two Higgs fields. CP violating complex particle masses varying 
with the Higgs field in the wall are essential. Since in the MSSM there is no 
indication of spontaneous CP violation around the critical temperature (contrary 
to the NMSSM) we have to rely on standard explicit CP violation. Using the 
WKB approximation for particles in the plasma we are led to Boltzmann trans- 
port equations for the difference of left-handed particles and their CP conjugates. 
This asymmetry is finally transformed into a baryon asymmetry by out of equi- 
librium sphaleron transitions in the symmetric phase. We solve the transport 
equations and find a baryon asymmetry depending mostly on the CP violating 
phases and the wall velocity. 
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1 Introduction 



The electroweak interactions violate baryon number. There is also a source of 
CP violation in the Standard Model (SM), and possibly even stronger ones in 
models beyond the SM, like supersymmetric theories. Thus it is very tempting 
to expect that these theories which are presently under stringent experimental 
tests allow to explain the creation of a baryon asymmetry in the early universe. 
However, in thermal equilibrium a baryon asymmetry (B + L ^ 0) is reduced by 
the electroweak interactions (conserving B — L). If there are no B — L violating in- 
teractions at the GUT scale creating B — L ^ 0, (as e.g. in the recently discussed 
"leptogenesis" models one needs a strong deviation from equilibrium to create 
a baryon asymmetry at the electroweak scale and fast freeze out in the subse- 
quent quasi equilibrium. Thus a strong first order phase transition with bubble 
nucleation is needed. However, combined analytical and lattice based work in the 
last years has made sure that such a phase transition is not present in the SM; 
actually there is no phase transition at all for Higgs masses beyond about the 
H^-boson mass 0. Thus, one has to inspect variants of the SM. Supersymmetric 
models are the most attractive ones. In particular the Minimal Supersymmetric 
Extension (MSSM) closest to the SM is investigated most intensively in exper- 
iments. With a light stop superpartner of the right handed top having a mass 
some GeV below the top mass, one can increase the 3 -term in the effective 
thermal Higgs potential leading to a strong first order phase transition [0-0]. 

In order to create a baryon asymmetry we first want to generate a chiral asym- 
metry in the hot phase in front of the bubble which is then transformed into 
a baryon asymmetry through hot electroweak sphalerons. For this we need CP 
violation in the Higgs field bubble wall or some CP violation in the hot plasma 
disturbed by the wall. As consequence, sources are induced in the transport 
equations. In the MSSM the chargino and the stop mass matrices contain ex- 
plicitly CP violating phases. The Higgs Lagrangian can obtain CP violating loop 
effects but we will also discuss in some detail the attractive possibility of a spon- 
taneously generated "transitional" CP violation just in the bubble wall around 



the critical temperature p|-|14f . Therefore general methods for the calculation of 
multi-scalar bubble walls are necessary. In Section Q we derive the relevant equa- 
tions of motion, including CP violation, and present in Section ^ a method to 
solve them in general. In Section |4] we present their applications to the considered 
supersymmetric models. 
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Further on, in Section [5] we will study the transport phenomena close to the 
bubble wall in the quasi classical particle "WKB" -approximation fl5HT8|. Thus 
we neglect soft boson classical field effects, which for example may play a role for 
the friction due a gauge field in the hot plasma [M], as well as off-shell memory 
effects to be discussed in the framework of Quantum-Boltzmann equations |2"0| - 
23]. The stationary velocity of the bubble wall is very important for the finally 
created baryon asymmetry. It is related to the friction term in the transport 



equations. This has been treated in ref. [23] and recently in [25] pointing to a 
very small velocity in the MSSM which, in general, is quite effective for producing 
an asymmetry. It turns out that it might even be too small for a maximum 
exploit [P^] . In ref. [27] a class of supersymmetric models with an additional 
gauge singlet has been worked out. These allow for a strong first order phase 
transition already at tree level. In our discussing of wall profiles and CP violation 
we will also include this interesting case. Concerning the transport equations in 
the WKB approximation we can return from the calculation in ref. [[28] including 
the singlet field to a critical discussion of the MSSM case which was already 
worked out in ref. |16| before and reconsidered recently |18 |. Finally we give 



numerical results for the baryon asymmetry in the given framework. 



2 Field Equations and CP Violation 

For our purposes it is sufficient to consider the one-loop corrected effective Higgs 
potential Vr with thermal masses. We use the standard notation for the MSSM 
[ p9|[ . In general, there are CP violating phases between the Higgs field moduli h\ 
and h 2 . Due to gauge invariance, the effective Higgs potential only depends on 
one combination 6 = 6\ + 6 2 of the phases, 



V T (H 1 ,H 2 ) = V T (h 1 ,h 2 , 



(2.1) 



Defining 



6 = 6,- 6 2 , (2.2) 
we can rewrite the kinetic terms for the Higgs bosons in the Lagrangian density 
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h\ + h\ 



hi -hj 



(2.3) 
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Since <%Vr = the Euler-Lagrange equation for 9 is 



{hi + h 2 2 )d»9 + (hj - hl)d»e = C 



(2.4) 



which minimizes the energy for Cq = 0. 

After elimination of 8^9 ( |2.4j ), the equations of motion for hi, h2 and 9 read 



d*d»h x + 



0. 



d^h 2 + 7T ^L_d,9d^9 + -i-V T 



(hl + h 2 ) 2 ^ 



h 2 h 2 

12 : d»9 



hl + h 2 2 



dh 2 



0. 



(2.5) 
(2.6) 
(2.7) 



Two limits are important. The first limit corresponds to a radially symmetric 
spatial solution of 



d 2 K 
dr 2 



2dh, 
r dr 



dV T {hj 
dhi 







(2.8) 



(without CP phase 9). The solution describes the initial state of nucleating bub- 
bles, the "critical bubble" . Eqs. ( |2.8|) are also to be solved for the determination 
of the tunneling probability |30|-p2]| . 

Second, constraining equations (|2.5| )-( pTT|) to a stationary wall ("domain wall") 
moving with velocity v w we are left with just one direction x = z — v w t at late time 
t perpendicular to the wall. This is the more important period for baryogenesis. 
With the assumption of stationarity and almost planar bubble walls the equations 
TBD-fO) reduce to 



h\h% 



(hi + hi) 2 
h 2 h\ 



9" 
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dh~! 
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{hi + hi) 2 dhr 



V t 
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h 2 h 2 

,L l' L 2 nl 

[hl + hf 



89 T 



0. 

0, 
0. 



(2.9) 
(2.10) 
(2.11) 



where the prime denotes d x . Vt is deformed by the plasma (moving in the wall 
frame) friction and has degenerate minima in the stationary case like Vt=t c - We 
roughly identify the deformed Vt with Vt=t c - 
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For the NMSSM singlet field S it is more convenient [|13|, Q to divide it up 
into real and imaginary components, 



S = n + ic, 



(2.12) 



implicitly introducing a phase variable Qg = arctan(n/c). The equations of mo- 
tion for n and c are of Klein-Gordon type 



d^n+^V T = 0, 
on 

d^c+^-V T = 0, 



(2.13) 



(2.14) 



where V T is the NMSSM potential |28[ 



In Sec. fO]we will see that especially the variations of the CP violating quanti- 
ties 9 and c in the bubble wall play a very important role in the generation of the 
baryon asymmetry. But first we will discuss more generally how these equations 
may be solved. 



3 How To Find Bubble Wall Profiles 

The first order phase transition is mediated by expanding bubbles. In order to 
compute the baryon asymmetry, we should follow the history of bubbles from 
the moment of nucleation, until the time when the broken phase fills the Uni- 
verse. After nucleation, there is in general a long period of stationary growth. 
Nucleation is characterized by time dependent solutions to the full equations of 
motion in a background of in general non-trivial profiles for temperature, veloc- 
ity and chemical potential. In the SM such solutions were investigated in [ j33f . 
We will not consider the full problem here but focus on the profiles of station- 
arily expanding almost planar bubble walls at the critical temperature T c which 
interpolate between two minima of the effective potential. This case is more in- 
teresting case for baryogenesis. In case of friction with the background plasma 
the corresponding temperature is the nucleation temperature T n p5 . 



For the bubble wall profiles in the MSSM without CP violation there exist 
numerical approaches to solve the problem of critical bubbles with two Higgs 
fields in |53fl and more general in |]35|] . In [jlTJ, the CP profile has been investigated 



in the background of a fixed Higgs profile. We will now discuss methods which 
enable us to find profiles in more general scenarios as shown in Sections |4.1| and 
4"72] as well as in [IT 
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The usual method for the SM case with only one Higgs field is to solve the 
corresponding single equation of type ( 2.13Q numerically by "turning around" 
the effective potential Vr and dealing x for a time t. Then the problem can be 
regarded as an initial value problem for the inverted potential — Vr- Usually an 
"overshooting-undershooting" procedure can be applied. This works well since 
there is only one direction in field space. Moreover, it can be implemented quite 
simply by the standard Runge-Kutta-method. 

The situation is completely different once there are additional directions in 
field space. Again one can consider the analogous mechanical problem with the 
turned around potential. Neglecting first derivatives which act like friction on 
the potential, the initial value problem is equivalent to a mass point rolling from 
the top of one hill (first minimum) along the ridge such that it comes to rest 
on the top of the second hill (second minimum). Small changes in the initial 
conditions lead to a completely different shape of the solution. In general it is 
not possible to know the initial conditions with sufficient accuracy to find the 
desired solution. Hence we have to device another method. Here we use the 
method of minimization of the functional of the squared equations of motion. 

A stationary problem for big bubbles can be reduced to a problem in one spatial 
coordinate x = z — v w t for a constant wall velocity v w . A constant velocity is 
caused by friction with the background plasma which ensures that the minima 
are degenerate |l,|2|,[g|32|,[36|. 

To solve Eqs. ( |2.9| )- ( p7TT| ) one has to find field configurations for which 

/+oo 
dx [Ef(x) + E 2 2 (x) + E 2 3 (x)] = 0. (3.1) 
-oo 

This is achieved by searching for the absolute minimum of S. Here Ei(x), E 2 (x) 
and E 3 (x) are the left hand sides of Eqs. (p.9|), ( |2.10| ), and ( |2.11| ), respectively. 
For the NMSSM the extension is straightforward 

/+oo 
dx [Ef{x) + El(x) + El(x) + El(x) + £ 5 2 (x)] , (3.2) 
-oo 

where E^(x) and E 5 (x) are the squares of stationary versions of the l.h.s. of 
Eqs. ( fn5D and (gTjj ): 

dVr 

E,{x)=n"-°-^- ) (3.3) 

dVr 

E 5 (x)=c"-^. (3.4) 
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This minimization method has also been successfully used in [K3] for the critical 
bubble. The approach of [57] is in principle also a minimization procedure and 
therefore the following considerations are also applicable. The method of |f?|.|3~5 



deals directly with the saddle point. Using the minimization method we have 
to solve a boundary value problem instead of an initial value problem. Thus we 
have to use an ansatz which fulfills the boundary conditions for each function for 
which we want to find the time development. 

Our procedure works in two steps. The first crucial step is to find an ansatz 
which is as close as possible to the exact solution. 



Step One: Ansatze In |6],|r^|34|,|38[] it was found that a kink is a good ansatz 



for the Higgs fields. We therefore use for N Higgs fields 



0f n * = f (l + tanh^ + x,))- / I ... \- (3.0) 



The Li and X{ are determined by minimizing ( |3.1|) or ( ]3.2|) which are built up 
from these configurations in the Eqs. of motion like ( |2.13|) -( 2.14|) . This first step 



is cheap in terms of computer time since it is only a iV-dimensional minimization 
and can be performed very fast. Minimizing with respect to a few parameters is 
a very successful first step since it reduces the value of ( |3.1| ) already significantly 
compared to a general function which only fulfills the boundary conditions. Hence 
the convergence is considerably improved. 

The shape of the tunneling trajectory depends strongly on the CP-odd Higgs 
mass parameter m^. Small values of give a larger mixing of the CP even 



Higgs fields and a sharper bending curve (see pg[|). 

Also the ridge between the the peaks turns out to be a rather fine ansatz. In 
general, the ridge is difficult to define but often it can be defined sufficiently in 
two scalar field dimensions. In a generic NMSSM case the ridge is even nearer to 
the solution than a superposition of two kinks. This is demonstrated in Fig. 0, 
where the solution and the ridge lie quite close to each other. 



ridge In two scalar dimensions an approximation to the ridge often can be 
defined as follows: Intersect the effective potential in planes perpendicular to the 
straight connection between the minima and take the maxima of the intersection 
curves in those intersection planes as points of the ridge. 

To use this in the NMSSM where we have basically 3 fields we define a new 
average field h = \Jh\ + h\ in the h\-h2 plane. This is a good approximation for 
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hi and /12 which are close to a straight line, similarly to the MSSM. They can 
subsequently be treated as one. Then, a good approximation to the ridge can be 
found by applying the above construction to the effective potential in the h — S 
space. In Sections ^4.2| we will discuss examples which violate CP. In these cases 
the ridge cannot be defined in the described simple manner. Fortunately these 
example cases again permit ansiitze utilizing the extended kink Q3.5| ). 



Step Two: Minimization and problems connected with it: With a pre- 
optimized ansatz we come to the second step, the high dimensional minimization. 
To represent ( |3.1| ) and ( |3.2| ) numerically we discretize them over a grid. We have 
to discretize the space variable x = z — v w t, the fields 0» (#_,-), the derivatives of the 
fields (first and second derivatives), and derivatives of explicitly given functions 
like the effective potential. These points have recently been discussed partially 
in |7|, |51| , P5[| . We want to discuss here some forthcoming details and computational 
problems. 

The fields are defined on a grid of points Xj along x. The first derivatives turn 
out to be appropriately discretized by f'{x) = (f(x + e) — f(x — e))/e and the 
second derivatives by f"(x) = (f(x + e)—2f(x) + f(x — e))/e 2 . Accuracy is of order 
e 2 where e is the grid spacing. The discretized field equations ( PTT| ) and ( |3.2| ) can 
be built out of these elements and the discretized fields. Derivatives including 
more points may even impair the result since they require more additions and 
subtractions which are numerically problematic. For this problem see also p9| , p|1 . 

Minimization parameters are the values <pi{xj) and 0(xj), respectively. With 
N fields defined at M grid points we have a, N x M-dimensional function to be 
minimized. 

In our applications the minimization itself is accomplished using two different 
methods, "Powell's method" |39[] and a "Pseudo-Langevin" method which will be 
discussed below. We also used the "downhill simplex" minimization method |]39] 
for comparison with which turned out to be too slow for practical use beyond 
a few computations. 

For the discussion of the practical usage of minimization we first discuss a 
common problem of minimizations including numerical derivatives. In general 
several undesired minima occur which can be categorized as follows: 

a) Some are real solutions of the equations of motion. A trivial example is 
9(x) = 0, which always solves (p. 5[)- (p77|) . 

b) Fake minima due to the numerical representation of the functional. This is 
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a common problem of discretized derivatives with finite differences. 

c) Minimization of ( |3.1| ) is achieved by solving 5S = 0, which is only a necessary 
condition for solution configurations. For a squared form, S = a 2 , also pseudo- 
solutions according to 5a = may exist in addition to the (desired) solution 



a = 0. Also combined effects may occur as discussed in p8 



Altogether this induces the problem that, starting from an ansatz which is in 
the vicinity of such an apparent solution, the algorithm might never converge to 
the desired exact solution. How can these problems be solved? 



One method is related to the method of "Simulated Annealing" |p8|,fn]|. There 
the trajectory through the high dimensional configuration space is given by a 
Langevin type equation, e.g. 

dqi dS[q] 

-aW=-^q- +7] ^ (3 - 6) 

where the qi are the degrees of freedom of the discretized action S[q] (Eq. (|3.f| ) 
or ( |3.2|) )and r]i(T n ) are random variables, typically Gaussian distributed. Every 
configuration is randomly varied to jump out of unwanted local minima. For a 
global progress in minimizing S a small suffer loss is permitted. An algorithmic 
temperature is defined which regulates the amplitudes of the fluctuations. This 
method is commonly used in lattice calculations (e.g. @). Unfortunately, the 
action 5* in our case is quite sensitive to randomly varying degrees of freedom 
through the derivatives of the configurations. A random change in the fields 
causes a strong change in the derivatives which increases S. The convergence is 
very poor and it turns out that this method does not solve our problem. 

Writing down the discretization of (|3.6| ) and omitting the fluctuation term we 
find the "Pseudo-Langevin" equation 

qi(r n+1 ) = qi(r n ) - At |^t . (3.7) 

Now the new time interval At of the new Hamiltonian system is a common 
integration step for all degrees of freedom. The motion promoted by Eq. Q3.7| ) is 
a motion without inertia through the configuration space. It converges towards 
a minimum of S. At is determined through the scale of dS[q]/dqi(n): 

A, ~ |M (3 ,) 
oqiin) 

In order to get a stable trajectory it must be sufficiently small. Now, a magnifi- 
cation of At after some time steps causes a "coherent" jump of the whole system 
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out of a local minimum. This is exactly the desired behavior of the "Simulated 
Annealing" method. But in that case the "decoherent" variation of the configu- 
ration shifts the system configuration to a completely unrelated point, permitting 
almost no effort at all in minimizing S. Thus, the system starts again from a 
worse configuration. It converges extremely slowly for a considerable number of 
degrees of freedom. However, the coherent annealing of our Pseudo-Langevin 
approach allows to "jump" over local hills without completely loosing properties 
of the configuration. It permits a larger time step to overcome small bumps and 
troughs, see Fig. |]. 



Figure V. Pseudo-Langevin method: Large time steps allow to overcome small 
bumps and troughs. Small time steps or simple minimum finding routines may 
lead into undesired local minima. 

This method is quite efficient and sometimes called "gradient descent" method. 
An interesting variant is also used in ]7|,[35| which directly deal with the saddle 
point configuration. 

The limit of our method is reached with derivatives ^|^y in (|3.7|) for the qi 
varying on strongly differing scales. Then the time step can not be computed 
appropriately for all qi simultaneously and the system runs over a minimum at 
one qi while another qj has not changed considerably at all. Such a behavior 
appears for large /i 3 -bumps in the effective potential. This may already be the 
case in the NMSSM for some parameter combinations. But a combination of all 
presented methods can also treat more complicated cases we cannot discuss here. 

What can be done to avoid some of the problems and to rate the quality of 
a minimum found? First it is important to have an ansatz as close as possible 




Starting point 
(ansatz) 



Small time steps Ax 



GlobaTminimum, S=0 
(solution) 
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to the desired solution to avoid reaching an unwanted local minimum after the 
time consuming second minimization step. One can also increase the number 
of intermediate steps in the procedure to improve the precision at the cost of 
computing time. Additionally we give an independent check to rate the results. 

Rating of solutions As long as we have energy conservation, one can check 
the quality of the results easily. The relation 

V[ qi ]/T[ qi ] = 1 (3.9) 
or V[qi] - T[ qi } = for all x (or r), (3.10) 

where T[gJ is the kinetic part of the action, must be fulfilled and one can check 
the deviation. With our Pseudo-Langevin approach the best solutions reach 
^[g]/T[g] — 1 with a precision of order 10 -3 . However, such a high precision is 
not necessary, even for sensible quantities like A/3. One can construct further 



identities that have to be fulfilled like an extended virial theorem as in |p7 . 

There is also a possibility for cases without energy conservation. Let us first 
discuss the problem of shooting along a given path in the Higgs field space because 
this method also gives a quality check for non energy conserving cases. Imagine 
that we have chosen a path for the bubble wall solution in field space. This does 
not mean that we have solved the equations of motion which would yield the 
x dependence of the fields. But for a given path we can use the overshooting- 
undershooting procedure to find the ^-dependence. The naive method to take 
the potential along the path as an effectively 1-dimensional potential gives wrong 
results even for cases with energy conservation. It can lead to a wall thickness 
which differs from the real result by a factor of two or more. The reason for 
this is that the differential equations only contain partial derivatives. In the 
naive approach the effective derivative (along the path) is a combinations of the 
original (partial) derivatives. But it can be extended to construct a quality check 
of solutions by looking at each field space direction separately. 

One necessary condition for reliable results is the shooting condition for the 
bounce solution: The "marble" rolling along an arbitrary path must come to 
rest at the zero minimum (symmetric phase). The field configuration q which 
defines the path coordinates is a combination of the scalar fields (hi, ti2, ■ ■ ■ )• A 
necessary condition for the solution configuration is 

dx^-^=0, (3.11) 
x , t oh k 
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which has to be fulfilled separately for all field directions k. Each configu- 
ration hk(x) can be projected to the chosen path q leading to configurations 
(jk(x), k = 1...N for N fields. Vice versa, assuming that q(x) is the exact 
solution to the (critical) bubble wall equations, all (jk(x) are then identical by 
construction. We can use this to construct a general check. The configuration 
q(x), obtained by shooting along a chosen path, may be projected to all com- 
ponents hk(x) in the field space which automatically fulfill ( |3.11| ) independently. 
The inverse functional x[h k ] can subsequently be used to construct a motion con- 
figuration q along the path for every k: qk{x[hk\). For the solution this leads to N 
identical configurations qk(x[hk\). But every pseudo solution different from the 
exact solution configuration leads to differences in these reconstructed pseudo- 
solutions which rates the quality of the pseudo solution. 
Concluding, the path has to be chosen that 

%end fly f^end fly 

dx—... = / dx— = ... = Vht (3.12) 



^ start 



dh\ J T t . dhi 

- 1 - *-* JL xf.n.rt. L 



(shooting condition) and simultaneously the (inverted) motion along the path 
must be the same for every field hi. This even works for critical bubbles where 
we have no energy conservation and the energy check does not work. 

The methods presented here can be used to investigate the dynamics of bubble 
wall expansion and baryogenesis in more detail. In Section |] we extend them 
to investigate CP violation for many parameter combinations. Our considera- 
tions and calculations were also essential to develop better perception of which 
approximations might be useful and appropriate to simplify calculations as done 
in Section EJ. More applications can be found in refs. |LT3],[25|,p7|,p5[ 



4 CP violation 



To calculate the actual baryon asymmetry CP violation as one of the Zakharov 
conditions becomes important. In WE, 



-IS | it was found that a sufficient 
baryon asymmetry requires CP violating phases as large as (9(10 _1 ). This might 
lead to contradictions with experimental constraints |46[] on the electric dipole 
moment of the neutron (EDM). It is still under discussion whether this contra- 
diction might be resolved by special cancellations. Large phases are also allowed 
if the 1st and 2nd generation squarks are heavy [fi7 . Regarding our results in 
the final Section |5] this becomes important again. For a more detailed discussion 
of this point see refs. MM and |BD] . 
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Figure 2: Solution and ansatze in the NMSSM. Left top: Fitted tanh(x/L + x)- 
ansatz and the ridge lying directly near the actual solution. Lower left: Shape of 
solution versus x compared to a kink ansatz. Right: 3- dimensional views of the 
effective potential with solution and kink ansatz. 

It is known, that in principle there may be spontaneous CP violation in the 
MSSM at zero temperature realized by a CP violating phase between the two 
Higgs doublets. While spontaneous CP violation at T = is forbidden for the 
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experimental allowed parameters ||, it was suggested that it might be realized 
at finite temperature quite naturally or especially only in the bubble wall 

between the symmetric and the broken phase []10, TT|, ^1] . 



An interesting scenario is therefore a temperature induced CP violation in 
the bubble wall during the first order electroweak phase transition which is not 
restricted by any experimental bounds. Hence it may be maximal and might 



sufficiently support the generation of the baryon asymmetry of the Universe flT7 , 
52|j . It will be referred to as transitional CP violation in contrast to CP violation in 
the broken minimum. In the literature p|- [l2| . ^0| .[5ll| the conditions for transitional 



CP violation were widely discussed. In Section [4.1| we will revisit this question by 



following the approach of ref. |13| in more detail. Since this approach is different 
to previous ones we compare the methods. 

Explicit CP violation is introduced by complex mass parameters in the effective 
potential. Independently from strong experimental bounds it might be interesting 
to study its effects on bubble walls and whether it suffices to contribute consider- 



ably to the baryon asymmetry. In Sec. |T2| we will analyze in detail the evolution 
of explicit CP violating phases in the bubble wall. 



4.1 Transitional CP Violation in the MSSM: Revisited 



In this section we elaborate on the systematic search of |H| for a parameter 
window for transitional CP violation in the MSSM. First we consider the case 
without explicit CP violating phases and check whether a CP conserving solution 
{9 = 0, 7r) is a local minimum of the action. Instability in the ^-direction requires 



ml{h 1 ,h 2 ) 



i PVTih^e) 



\h1h2 



de 2 



< 0. 



(4.1) 



0=0 



h<2) is a measure for the instability at a given point in the wall. Vr is the 
finite temperature effective potential including stops, charginos and neutralinos. 
In our convention hi can have either sign which permits to consider only 6 = 0. 
Condition (|4.1|) agrees with the constraint of Lee |53| on which |9|, pi, fall are 
based on. But Eq. (|4.1| ) is true more generally along the whole bubble wall. 

We are interested in regions where ( |4.1|) is fulfilled. In our search we first 
neglect the strong experimental bounds on the Higgs mass since it might still 
be interesting to find a shape of such a transitional CP violation. The region of 
generically allowed parameters is an 8-dimensional "hypercube" stretched by the 
tree level parameters itia and tan/3, the stop mass parameters tuq, mu, trilinear 
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couplings At, and fi as well as the gaugino mass parameters M\ and M 2 in the 
neutralino and chargino mass eigenvalues. For each chosen parameter set we must 
determine the critical bubble wall profile h\(z), h 2 (z) including the 9(z) field and 
scan the wall for points which fulfill the instability condition ( |4.1| ). 

The methods for solving the (critical) bubble equations with more than one 
scalar field including 9 were discussed in Section |] and in . However, all these 
methods are very time consuming and need sophisticated algorithms. There- 
fore we proceed in two steps. First we will restrict the possible parameter space 
in order to obtain a more promising parameter set permitting transitional CP 
violation. Hence we will investigate the averaged dependence on the parame- 
ters. Second, with this preparation we will search for a definite parameter set 
permitting the desired CP violation in the bubble wall. 

A first and useful simplification is to refrain from determining the wall profile 
very accurately for each parameter set. Indeed, in the MSSM the solution to the 
equations of motion is quite close to a kink ansatz h(z) ~ l/2(l+tanh(z/L)), and 
therefore we do not need the full solution to the equations of motion. The solution 
can be approximated by the ridge of the negative potential —V e ff. Actually, it is 
even sufficient to approximate the solution in field space h 12 by the straight line 
between the minima. This approach is correct as long as there is only small CP 
violation which does not have a too strong back-reaction to the wall shape as could 
happen for large phases (see ref. f38| ). If we found a candidate for transitional CP 
violation with a large phase, we would have to solve the whole set of equations. 
But for the moment we are only searching for promising parameter regions. 

In order to perform a complete search for a parameter window with transitional 
CP violation in the bubble wall it is important to consider the dependences on 
the various parameters. At tree level we find 

mj 2 = ~m 2 A sm2{3. (4.2) 

Hence, the minimum in 9 direction is at 9 = 0. Thus, in order to obtain sponta- 
neous CP violation the 1-loop contributions have to overcome the tree level term. 
Obviously large tan/3 and small tua are preferred. The dependence on other pa- 
rameters is not quite obvious. Potentially it is possible that m 2 is negative for 
very small (negative) mjj « (2itT) 2 « m 2 Q |T2|,p|. 

In an analytical approach where we plot m 2 versus various parameters we 
always fix the remaining parameters. We find that it dominantly depends on 
rriA and tan/?. All other parameters are effectively 1-loop corrections and give a 



14 



correspondingly smaller effect. The strongest 1-loop dependence is given by A t 
and \x (see Fig. |3|). We also find here that large values with opposite signs give 
the smallest m|. Later we will see that the latter effect appears only to be correct 
in these preliminary analytical investigations. 




800 800 



Figure 3: Dominant 1-loop effect: in dependence of A t and \i including stops 
and charginos. The tree level contribution causes a shift in the mean 777J. Large 
values with equal signs are preferred to give a negative contribution to rnj. 

An interesting result we obtain by this method is that mj(/ti,/i2) tends to 
become larger while approaching the broken phase within the wall. The opposite 
behavior was noticed in [[II]]. An increasing m\ is undesired since large phases 
in the broken minimum are restricted by the experiments. Also small phases are 
excluded by experimental bounds on the Higgs masses which are connected with 
small CP violation ||. So even if a transitional CP phase occurs analytically it 
might be excluded by experiments. 

The other parameters give minor effects which nevertheless can contribute. 
Large M 2 and [i are favored in order to get small mf. As already mentioned, in 
this kind of analysis we always fix most of the parameter and vary only some of 
them. But we are interested in the whole parameter space. We will see that the 
shown dependences must be treated with caution since they may change using 
more general investigation methods. 
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In the next stage we want to investigate the complete parameter space in a wide 
range. We therefore choose a different strategy. We will analyses correlations in 
average. At the first stage, we do not solve for hi,h,2,T c , but rather take them 
as free parameters in the ranges h\/T = —2. ..2 and hijT = 0...2, T = 80. ..120 
GeV. The zero temperature parameters are varied in the wide ranges 

tan/5 = 2.. .20, m A = 0...400 GeV, 

m v = -50.. .800 GeV, m Q = 50.. .800 GeV, (4.3) 
fx,A t ,M 1 ,M 2 = -800.. .800 GeV. 

Here a negative m\j means in fact a negative right-handed stop mass parameter, 
— \m\j\. We have also studied separately the (dangerous f7|) region where the 
transition is very strong |3|,|5| corresponding to mu ~ —70... — 50 GeV. 

Note that since we do not solve for the equations of motion at this stage but 
allow for hi = ±|/ii|, we have to divide in (^4.1| ) by hih2 instead of \hih2\~- this 
leads in general to positive values due to the tree-level form of the potential. 
A signal of a potentially promising region is then a small absolute value of the 
result, since this means that we are close to a point where dgVr{hi, h 2 , 9) crosses 
zero. 

In Fig. H the relative number of m| versus its occurring values demonstrates 
which values are the most likely ones. There is a strong peak at a certain positive 
value and the probability decreases rapidly for zero or negative values of mf . But 
nevertheless there are some negative m\. It has to be checked, whether they are 
relevant for the bubble wall. 

We found that spontaneous CP violation is statistically independent of the 
actual critical temperature. For checking the individual cases it may nevertheless 
play a role [T^,|TB . 



In the analytical result we observed that large /i and A t with equal signs are 
preferred. This can be checked in a two dimensional correlation plot in Fig. || 
Both combinations are correlated with large positive m|, hence with in average 
more stable regions. The best region is a cross-like region where either A t or fx 
are small. This could not be expected from the analytical approach which showed 
a preference of large values of A t and \i. 

The strongest correlations come from A t and /i. The scale of their variation 
is larger by a factor of about 3 than the scales of the 1-loop mass parameter 
dependences with averaged 1-loop parameters. 

Unfortunately small itia and large tan/? are strongly preferred to obtain an 
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Figure 4: Relative number versus value of m\. Small positive values are preferred. 
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Figure 5: Correlation of m\ and the Higgs mass tuh- Most favorable values are 
correlated with unphysical small Higgs values. 



unstable CP direction. This is in contradiction to the requirements of a strong 
first order phase transition ||. Moreover CP odd Higgs masses up to itla ~ 
100 GeV are experimentally excluded. 



17 



In [|J it was argued that spontaneous CP violation at T = is excluded be- 
cause of the experimental exclusion of the corresponding small Higgs masses. Ac- 
cordingly, for transitional CP violation physical Higgs masses lower than around 
40GeV are fairly probable, see Fig. ^. Here the average of m\ decreases rapidly. 




is the cross like region around vanishing A t and /x. 

There are also minor effects from the remaining parameters. The stop param- 
eters mu and m§ are preferred to be small (Fig. 0a). The effect is smaller than 
the main 1-loop contribution from A t and \i. Mi and M 2 are almost completely 
uncorrelated. There are smaller correlations with A t or /i (Figs. 0b, c) which are 
folded into a vanishing correlation of m\ with Mi, Mi (Fig. 0d). The effects 
on mu, rriQ, Mi, and M 2 do not strongly promote a small or negative m\, but 
nevertheless they can be taken as hints. 

In a second stage, we study the relevant parameter region in more detail. We 
start from the full parameter range ( |4.3| ) and restrict the search to the most fa- 
vorable region. First, we determine the critical temperature. Then the equations 
of motion are solved for (hi, hi). By comparing with the numerical exact results 
in several cases we find that a sufficient accuracy can be obtained by searching 
the "ridge" as an approximation, see Section [3]. Finally, we look for the mini- 
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Figure 7: Correlations of the 1-loop parameters in average. The measure 777J for 
instability versus the average of the other parameters is shown. The strongest 
correlation appears for the mass parameter /i. M\ and M<i show less average 
dependences. Small mu are preferred, which coincides with the requirements of 
a strong phase transition whereas also small rriQ are preferred, which points into 
the wrong direction. 

mum in 9 at fixed (hi, h 2 ) which is a fast and reliable approximation to the real 
solution as long as 9 is small. We make some further restrictions: we exclude 
unphysical negative mass parameters as m~ < and discarded extremely weak 
phase transitions v /T <C 1. We also excluded cases leading to T = spontaneous 
CP violation in the broken phase since this requires very small values of vtia ||- 
While in the statistical analysis we included 10 9 parameter sets, we investigated 
~ 2 x 10 6 individual configurations in more detail. Nevertheless, we could not 
find any case with the desired property of temperature induced transitional CP 
violation within the corresponding bubble wall in the MSSM. Furthermore, in 



19 



11|] a special point around mfj rs was considered. Since thermal mass were 



neglected, this corresponds to the point were m\j + #T 2 « 0. This point is in 
the vicinity of charge and color breaking minima. Without expanding the 1-loop 
contribution in V\/v2, we cannot reproduce the instable behavior proposed in [ |TT 



Moreover we found that in the MSSM the dominating effects which support 
small m\ are in contradiction to strong phase transitions and baryogenesis re- 
quirements. 

4.2 Explicit CP Violation 

We investigated also explicit CP violation in the wall. Although there are strong 
experimental constraints on the magnitude of the corresponding phases it might 
be interesting to see how strong the variation of 9 along the phase boundary can 
be, and if it contributes considerably to the baryon asymmetry. 

There are several possibilities to include explicit CP violation in the effective 
potential. In principle, several mass parameters might be complex. At tree level 
this is m\ 2 , and at 1-loop level we have more possibilities through the masses of 
stop, chargino, and neutralino. In the stop, chargino, and neutralino sector /i, A t 
and the soft breaking mass terms can carry phases: 

9 mi2 ^ m\ 2 = fiB, 9 A t <r+ A t , Optr* n, 

6 Ml <-> M x , 6 M2 <-> M 2 . (4.4) 

Altogether we have five phases. 9 mi2 can be absorbed in the Higgs fields and then 
arises in the mass matrices. The phases in the stop off diagonal entries can be 
absorbed into a complex A t — A t — fi cot /3. We define the remaining three phases 
as 6 A = 9 mi2 + 9 A t + 9^, 9 C = 9 mi2 +9^ + 9 M2 , and 9 M = 9 mi2 +9^ + 9 Ml where 
9a appears in the stop matrix, 9c in the chargino matrix, and together with 9m\ 
in the neutralino matrix. In the chargino and stop eigenvalues we then have to 
replace cos 9 — > cos(6 ) + 9a) and cos^ — > cos(6 l + 9c)- The replacement in the 
neutralino eigenvalues is more complicated and will not be given here. 

Again we investigate the variation of 9 along the wall. As in the spontaneous 
case, we search for minima along the 9 direction. For baryogenesis, we need a 
strong variation of the phase d x 9(x) along the wall. But even for maximal angles 
9a = 7r/2, 9q = 7r/2, and 9m 1 = we find a strongly suppressed CP phase 
in the broken Higgs phase. For itla < 90Ge\^ 9(x) is of order 10~ 3 . . . 10~ 2 and 
varies only moderately within the wall as shown in Fig. Especially there is no 
kind of peak as proposed earlier . Only for experimentally excluded values of 
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171 a ^> lOGeV we do obtain phases up to order unity. For explicit phases of order 
1CT 1 the dynamical phase is typically of order 10~ 4 . We again realize that the 
size of the variation is dominated by the tree level parameters and tan/3. 

From a technical point of view the evolution of the phase 9(x) into the sym- 
metric phase is of interest. It is more difficult to determine (and less meaningful). 
The solutions of Fig. |]have been obtained by using a tan/i-ansatz also for 6(x), 
which turned out to compare very well with the numerical precision solution in 
the middle of the wall, where the full equations of motion can be solved more 
easily. 




Figure 8: Dependence of 9 and 9' on itla near physically realistic cases: The 
derivative strongly depends on and tan/3. Solid: = 80GeV, tan/3 = 2.0, 
dashed: uia = 120GeV ', tan/? = 2.0, dot-dashed: uia = !20GeV, tan/? = 3.0. (x 
and 9' in units of GeV^ 1 .) 



4.3 CP Violation in the NMSSM 

In the NMSSM we have many possibilities and a much larger parameter space 
[ p 7 ^ . Besides the Higgs fields hi and h 2 we have a singlet field S. Permitting 
complex fields, we can write the singlet field S = n + ic implicitly introducing 
a phase 9 S = arctan(c/n). Indeed, it turns out that transitional CP violation 
can occur in the NMSSM quite naturally. In |13j an example is shown where in 
the SU(2) x L r (l)-symmetric phase the singlet carries a vev in both components: 
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(n,c) symm . = (50AGeV,99GeV) and (n,c) brokcn = (-122GeV,0GeV). Moreover 
the standard Higgs phase varies between zero and a value of order 9 sym ~ 1/20 
(in this special example) in the symmetric phase and zero in the broken phase. 
Both phases, 9 and 9 S start with non-vanishing values in the symmetric phase 
and vanish in the broken phase. 

We used the values k = 0.4, A = 0.05, /i = 212.6GeV, m Q = 278.1GeV, 
mu = 209.0GeV, B = 83.9, r = 533.2, tan/5 = -5.0, A t = -219.8GeV, A k = 
50GeV, A\ = -150GeV, M x = 41.3GeV, M 2 = 82.2GeV The critical 

temperature is T c = 102.2GeV. 

This demonstrates the general possibility of transitional CP violation in bubble 
walls of the NMSSM phase transition. We have a strong phase transition in the 
example case. This shows the principle viability of baryogenesis in the NMSSM. 
More detailed calculations can be found in I fZHf . 



5 The Baryon asymmetry in the MSSM in the 
semi-classical limit 

5.1 Generalities 

As just discussed there can be explicit CP violation in the MSSM and also tran- 
sitional CP violation in the NMSSM. The expanding bubble walls of a first or- 
der phase transition sweeping through the hot plasma generate temporary non- 
equilibrium. In the hot plasma we have sizeable baryon number violation by 
the weak sphaleron which rapidly stops in the equilibrium Higgs phase. Thus 
baryogenesis in principle is possible. However, we have to describe a concrete 
procedure since these are only necessary conditions. Technically one has to study 
transport equations in the plasma for particles/fields with varying mass due to 
the varying Higgs fields in the bubble wall. For thermal momenta p ~ T and a 
"thick" wall (L w ^> 1/T) the semi-classical approximation should be appropriate 
and we should end up with Boltzmann transport equations. In principle one 
has to start with the Baym-Kadanoff equations of the real time non-equilibrium 
formalism and end up with classical particle equations (plus some corrections 
as memory effects). In a different corner for soft bosonic fields a classical field 
description is expected. Here we only will deal with the particle limit and discuss 
the WKB approximation in accordance with refs. [16[ 18[| . Since we are mainly 



interested in the baryon asymmetry we will first calculate the transport for an 
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asymmetry between left-handed particles and their CP conjugates and then con- 
sider baryogenesis via (out of equilibrium) sphaleron transitions in the symmetric 
phase. The transport equations for particles alone are interesting if the friction 
of the bubble wall and the resulting stationary wall velocity is calculated. Here 
CP violation can be neglected. As it will turn out that the wall velocity strongly 
influences baryogenesis this is an important point, not being dealt with here, but 
we refer to refs. |T9|,p4|-|2l 



5.2 WKB approximation and dispersion relations 

In the following we summarize the derivation of the chargino dispersion relations 
in the presence of a slowly varying Higgs field background. For the details of the 
calculation we refer the reader to ref. where also the application to squarks 
was worked out. Indeed the contribution of squarks turned out to be sub-leading. 

To obtain the dispersion relation we solve the Dirac equation for the charginos 
in the WKB approximation 

ia^ -Aft J k 



0, (5.1) 



where 



r0\* 



Af-[ M2 92{H °^ ) (5 2) 



and ^ l = (W + ,h.2) T , = (W~ , hi) T . Since spontaneous CP violation does 
not occur in the MSSM, we have to rely on explicit CP violation by the M 2 
and \x parameters. At this level all interactions between the particle and the 
plasma are neglected. The scattering effects will be accounted for in Section 5.3 
when Boltzmann equations are written down to describe the local phase space 
distributions. 

During the passage of the bubble wall the chargino mass matrix becomes space- 
time dependent. We neglect the curvature of the wall and boost to the rest frame 
of the bubble wall. Then the mass matrix only depends on the position relative 
to the wall, i.e. M = M(z). 

Exploiting conservation of energy and boosting to the Lorentz frame where 
the momentum parallel to the wall vanishes, we take the ansatz \1/ = e~ lEt ^(z). 
Eq. fl5J] ) separates into two equations |L7] 



id z t ± = ±Q{z)Ct, Q(z) = [ M E {z) M ^ Z) ) (5.3) 
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where £± are the S z = ±~ components of £. To solve Eq. ( |5.3| ) one brings the 
^-dependent matrix Q(z) = D(z)Q£ ) (z)Z)(z)~ 1 into a diagonal form, where p8| 



VE 2 



2 



The SU(2) matrices 



U \ f coshX e^sinhX 
V [ e^sinhX coshX 



jj _ [ cos a e n sin a \ _ / cos 6 e 1(5 sin 6 \ 

' — e n sin a cos a ) ' \ — e" 5 sin 6 cos 6 ' 



are determined by M = VM^U^, where Mp = diag(mie* 6 ' 1 , r^e* 92 ). The entries 
of X which is a diagonal 2x2 matrix in flavor space, obey tanh2X 12 = ra\ ijE. 
In the local WKB basis £± = D~ l ^± the Dirac equation ( p.3|) takes the form 

ihd z £ ± = (±Q D - D- l ihd z D)l ± (5.6) 

which still is an exact equation. For mixing Dirac fermions like the charginos 
the transformation to the WKB basis consists of a rotation in flavor space (first 
factor in D) and a rotation to the helicity basis (second factor in D). In case of a 



single Dirac fermion, e.g. the top quark, there is only the helicity rotation [17| ,g8], 
while for scalar particles, e.g. the stops, only the flavor rotation is present [p§[| . 

In general, the correction term D~ l ihd z D caused by the position dependent 
field redefinition is not of diagonal form. The components of £± are still coupled. 
However, the off-diagonal part is proportional to d z D ~ D/L w . Typical momenta 
of the particles in the plasma are of the order of the temperature T, which is 
much larger than 1 / L w for the bubbles under consideration. We therefore expand 
Eq. ( |5.6|) in powers of d z or more precisely in powers of h (WKB approximation) 
that we already reintroduced for that reason. 

To order (h)° we can neglect the D~ l hid z D contribution. Thus the components 
of £± decouple in Eq. Q5.6Q. Inserting the WKB ansatz for the fermion field 



into ( |5.6|) , we obtain the dispersion relations 

p z (E) = sgn(p z )V 'E 2 - m 2 . (5.8) 
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Obviously, in the classical limit the ^-dependence completely disappears, demon- 
strating that CP violation is indeed a quantum-mechanical phenomenon. 

To solve the Dirac equation ( |5.6|) to order h we have to take into account the 
D~ l Tiid z D term which reintroduces a coupling between the two components of 
The dispersion relations p z {E) are obtained from the eigenvalues of the matrix 
±Qd — D~Hhd z D. Since to order % the off-diagonal terms do not contribute we 



are left with 16.28 



Lj : p z = sgn(p 2 
Li : p z = sgn(p 2 
Rj : p z = sga(p z 
-Rr : p z = sgn(p 2 
where 7=1,2 and Q\ = d z 9i, etc. and 
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— rrij 


E 2 
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— rrij 


E 2 
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— rrij 


E 2 
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— rrij 



% + 5' sin 2 b) sinh 2 X/ + 7' sin 2 a cosh 2 X/, 
9'j + 5' sin 2 6) sinh 2 X/ — 7' sin 2 a cosh 2 X/, 



E 2 — m 2 - + (ff t - 7' sin 2 a) sinh 2 X 1 + 5' sin 2 6 cosh 2 X T , 



sinh 2 X 7 = (E - VE 2 -m 2 )/2v r E 1 



m 2 . 



(5.10) 



In the symmetric phase L 2 and R 2 evolve to the left-handed Higgsinos states 
h 2 and hj, respectively. The flavor transformations U and V are related to the 



parameters of the chargino mass matrix ( p72 



sin 2 a = 2|A| 2 /A(A + A) with 

a = s 2 ((M 2 /f 2 °r + 



A 

A : 



\M, 



1 2 



(5.11; 



(A 2 + 4|A| 2 ) 1/a 
and 7 = arg(A). This gives 

i sin 2 a = 2Im{A*A') /A (A + A) 

and there are similar relations for sin 2 b, 5 and 5' sin 2 6 exchanging a and 6, 
and H 2 , 7 and —5. The mass eigenvalues (in non-symmetric notation) read 



[M D )u 
[M D ] 22 



, , cos a 

M 2 ^ 

coso 

cos a 



cos 6 



(5.12) 



cos 6 ^ v LJ cos 6 

Notice that the variation of mi, which is encoded in d z Xi, d z a and d z b drops 
in the dispersion relations. The CP violating part of the dispersion relation is 
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proportional to the derivative of the phases 6j, 7 and 5. Thus only varying phases 
contribute to CP violation in the semi-classical limit. Furthermore, CP violation 
is proportional to derivatives which guarantees that its effect is turned off far away 
from the bubble wall. Because of the different dispersion relations, left- and right- 
handed particles feel a different (semi-classical) force in their interaction with the 
wall. For the anti-particles, L and R, the CP violating part comes with the 
opposite sign. 

As pointed out very recently in ref. [IS] one should use the kinetic momentum 
Pkin = mVgroup = rn^ instead of the canonical momentum p (we used up to now) 
in the classical limit leading to Boltzmann equations, quite in the spirit of the 
correspondence principle of basic quantum mechanics. Calculating the (inverse of 
the) group velocity from ( |5.9| ) ^-independent terms drop out, i.e. cosh 2 Xj can be 
substituted by sinh 2 Xj. The kinetic momenta beyond the zeroth order (5.8) then 



contain correction terms ±(#j + <5'sin 2 6 — 7' sin 2 a)m 2 / (2E(E 2 — m 2 ) 1 / 2 ). Thus 
indeed kinetic and canonical momentum are not equal because of CP violating 
effects. 

The dispersion relations can be inverted. To leading order in the derivatives 
the CP violating part of the dispersion relation for the eigenstate L 2 which in the 
symmetric phase corresponds to h 2 is 

AE = -sgn(p z )(6' 2 + 5 , sin 2 6-7 , sin 2 a)m 2 /2(^ in + m 2 ). (5.13) 

(For pkm ^> Tn this is twice the result one would obtain with the canonical mo- 
mentum after substituting cosh 2 by sinh 2 in ( f5.9| ).) Most importantly AE is now 
totally symmetric under the exchange of Hi and H 2 . This will destroy the most 
prominent source term ~ H\H' 2 — H[H 2 of older work. 

The phases 7 and 5 only vary due to a change in the Higgs vev ratio tan/? or 
because of transitional CP violation in the bubble wall. The first contribution is 
suppressed, since the variation of (3 is at most ~ 1CT 2 for realistic Higgs masses 
[B4|, B3], while transitional CP violation most probably does not occur at all in 



the MSSM as discussed in Section 4.1 ||13|| . On the other hand, the contribution 



to the chargino dispersion relations stemming from the variation of the complex 
phases in Mo only requires explicit CP violating phases in fi or M 2 . Eq. ( |5.12j ) 
demonstrates that even though the phases in the two terms entering [M D ] 112 2 
are position independent, their contribution to the resulting phase varies due to 
the change in the (real) Higgs vevs. 
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5.3 Diffusion equations and the baryon asymmetry 

In this Section we study the coupled differential equations that describe particle 
interactions and transport during the phase transition. We treat the plasma as 
consisting of quasi-classical particles with definite canonical position and momen- 
tum. The phase space distributions fi(x,p,t) of the particles evolve according to 
the classical Boltzmann equation 

d t fi = {d t +'x-d s +'p- d p )fi = d[f\. (5.14) 

The time derivatives of position and momentum obey the Hamilton equations 
x = dpE(x,p) and p = —d$E(x,p), where E(x,p) are the dispersion relations 
derived in the previous section. This treatment is an approximation to quantum 
Boltzmann equations which have to be discussed in principle. This picture is 
justified for thick walls (p ^> 1/L W ) if it predicts a sizable effect, not dominated 
by non-leading terms in the derivative expansion |20[ . 



The Boltzmann equation can in principle be solved numerically. However, to 
make it analytically tractable we use the fluid-type truncation fl] 



fi&P,t) = em _ ViP ^ i)±1 (5.15) 

for the phase space densities of fermions (+) and bosons (-) in the rest frame 
of the plasma. Here and /ij denote the velocity perturbations and chemical 
potentials for each fluid, respectively. We also split E$ into a dominant part 
Eoi = \J p 2 + Tnf and a perturbation AEi ~ d z 6 which is related to CP violation. 
The chemical potentials are the central quantities that finally will determine 
the baryon asymmetry. The velocity perturbation on the other hand, is only 
introduced to allow the particles to move in response to the force, giving rise to 
chemical potential perturbations. 

We are looking for a "stationary" solution of the Boltzmann equation, because 
at late times the wall moves with constant velocity v w . This means that any 
explicit time dependence enters in the combination z = z — v w t. Inserting the 
fluid ansatz into the Boltzmann equation ( |5.14j ), eliminating i>j and taking the 
difference between particles and anti-particles, we obtain to linear order in the 



perturbations AEi and and to leading order in the wall velocity JT6 

- K^Difx" + VuHi) + ^Tp^fij = Si, 

p j 

s - ^ wk iP ^' (5 - 16) 
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where primes denote d 2 . The diffusion constants read Di = Ki(p 2 z / E )l/ (p 2 z T^), 
and T e , T d are the rates for elastic and inelastic processes, respectively. The 
averages are carried out according to 

n _ W4(-) _ j#pm (517) 

Here f'± = df±/dEo = — /3e /3jE °/(e^ £ ' ± l) 2 is the derivative of the unperturbed 
Fermi-Dirac or Bose-Einstein distribution. The statistical factor k is 1 for mass- 
less fermions, 2 for massless bosons and exponentially small for particles much 
heavier than T. The subscript "0" denotes averaging with the massless, unper- 
turbed Fermi-Dirac distribution. 

The CP violating source term Si in (|5.16| ), which is due to the semi-classical 



force, is proportional to the diffusion constant. This is because particles must 
move in order to build up perturbations. In Si the thermal averages over the CP 
violating energy perturbations AE are performed using the massive distribution 
functions in order to account for Boltzmann suppression of heavy particles. 

We obtain the transport equations of the MSSM from Eq. ( |5.16| ) by specifying 
the relevant particle species and interactions in the hot plasma. The network of 
equations can be simplified considerably by using conservation laws and neglect- 
ing interactions that are slow compared to the relevant time scale, which requires 

DK « r-i. 

In a first step we neglect the weak sphaleron interaction with rate T ws , which 
will be included at the end of the calculation. In the following we will there- 
fore assume baryon and lepton number conservation. The neglect of the weak 
sphalerons allows us to completely forget about leptons in our transport equations 
and compute only the quark and Higgs densities. 

We assume the supergauge interactions to be in equilibrium. The chemical po- 
tential of any particle is then equal to that of its superpartner, and it is convenient 
to define the chemical potentials (iu = (/i 1i c + /i flc )/2, (Iq 1 = (/i u + /i<i + + j)/4, 
A* Hi = (/%o + fJ-H- + A^o + /i^-)/4, etc. Furthermore, we take into account the 
following interactions 

(T y + r yA )(fi H2 + fi Qa + fir), T yfl (fi Hl - fiQ 3 - fpp), 
r ss (2/ig 3 + 2/iq 2 + 2/i Ql + fir + fiB + A*c + + Vu + A*d), 
^/(/Affi + fJ'H 2 )i r m (/ig 3 + (A T ), T Hl (i Hl , Yh 2 (i H2 . (5.18) 
The rates in the first line are related to the Lagrangian 
Ant = ytt c qsH 2 + y t t c q 3 h 2 + y t t c q 3 h 2 - ytfd^^Hi + y t A t t c q 3 H 2 + h.c. (5.19) 
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T ss denotes the strong sphaleron rate. Thf is due to Higgsino helicity flips induced 
by the fxhih 2 term. T Hl 2 and T m correspond to Higgs and axial top number 
violating processes, present only in the phase boundary and the broken phase. 

If the system is near thermal equilibrium, number densities and chemical po- 
tentials are related by 

m = IkifaT 2 (5.20) 
o 



where ki is the appropriate sum over statistical factors k introduced in (15.171) , 

e -g- k Ql = N C (K U +K d +Ku+Kj), k v = N c (k u o+Kuo), k Hl = (K H o+K H -+K- h o+K- h -), 

etc. N c = 3 denotes the number of colors. In the massless limit one obtains 
^<2i,2,3 = 18 > k v = k D = ... = k T = 9, k Hl2 = 6. 

Using baryon number conservation and neglecting the small Yukawa couplings 
of the first and second family quarks the strong sphaleron rate reads 



T ss (2fiQ 3 H h ix d) = T s . 



2 + 9^)^(1-9!)^ 



(5.21) 



To arrive at this expression we made the assumption that all the squark partners 
of the light quarks are degenerate in mass. Assuming equilibrium for the strong 
sphalerons we obtain 

2k B + 9k Qs 

^ = 9k T - k B (5 - 22) 

The reduced set of diffusion equations for the relevant particle species then finally 
reads 



- AV q fi Q3 + (T y + Y yA )[ix H2 + Bii Q3 ] - T vlt \fi Hl - Bii Qa ] + BT m fx Q3 = 

— kaJ^hlXHx + ^yulP-Hi — BfXQ 3 ] + ThfiUHi + l i H 2 ) + ^HifJ-Hi = ^Hi 

-k H2 V h ix H2 + (T y + T yA )\p H2 + Bix Q . 3 ] + T hf (ix Hl + IXh 2 ) + ^h 2 ^h 2 = Sh 2 (5.23) 
where 

9k T k Q3 + 9k B k Q3 + Ak B k T 



A 



9kx — k B 



B = kB + 9kT + 9kQ > (5.24) 
9kr — k B 

andP, = Di& + v w £. 

We keep the rates related to the top Yukawa interactions finite. If these in- 
teraction are in equilibrium, the resulting diffusion equations are sourced only 
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by the combination Shx — Sh 2 , because of the constraint + /xh 2 = 0. As a 
result the dominant contribution to the chargino source terms stemming from 
the ^-dependent part in the dispersion relations Q5.9Q cancels, because the cor- 
responding terms for h{ and h% are exactly of the same size. This is not true 
for the 7', 5' contributions which arise from the "flavor" transformations U and 
V. However, they are suppressed by the small variation of tan/3 in the bubble 
wall §,0,0. 



The full diffusion equations ( 5.23|) have already been studied in ref. [P8], with 



source terms corresponding to the helicity part of the dispersion relations (|5l 
In the following we also include the flavor part of the source term in the analysis. 

Before solving the network of diffusion equations, we turn to baryon number 
generation by the weak sphaleron processes. The evolution of the baryon number 
density n B is governed by 

-V q n B + 3Q{z)T ws {T 2 fi BL - an B ) = 0, (5.25) 

where we have assumed identical diffusion constants for all quarks and squarks, 
and neglected contributions of leptons. The position dependence of the weak 
sphaleron rate is modeled by a step function 0(z): anomalous baryon number 
violation is unsuppressed in the symmetric phase (z > 0) and suddenly switched 
off in the broken phase (z < 0). Baryon number generation is sourced by the 
chemical potential of left-handed quarks 



k Q ~ + 2k T (2k B 2k 



9kr — kn \ k 



+ 



T — l%B \ tQi tvQ 2 



HQs- (5-26) 



The second term in Eq. ( |5.25 ) describes damping of the baryon asymmetry by 



weak sphalerons in the symmetric phase. The parameter a depends on the degrees 
of freedom present in the hot plasma. Taking only the right-handed stop to be 
light gives a = 48/7 0. 

From Eq. ( |5.25|) one can easily obtain the baryon to entropy ratio in the broken 
phase 

f]B = — ^f^' r d ^B L (z)e-* (5.27) 







s 2'K 2 g Jf v w T 

where we have taken the entropy density s = (27r 2 g,,,/45)T 3 and v = 3aT ws /(2v. 



lq]. g* ~ 126 is the effective number of degrees of freedom at the phase transi- 



tion temperature. Eq. (|5.27| ) shows that the integral over the left-handed quark 
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number, til oc in the symmetric phase determines the final baryon asymme- 
try 

We now return to eqs. (|5.23|) in order to compute \xb l - These linear second order 
differential equations can be solved by finding the appropriate Greens function. 
In contrast to ref. we do not make any approximations in computing this 
Greens function. We keep the discussion general and consider the following set 
of N coupled diffusion equations 



/ -k n Vu + r 
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—ki^ViN + Tin \ 



\ 



—k N {D N i + r 



NX 



( Ml \ 



-k NN V NN + T NN J y fi N J 



( s 1 \ 



\S N ) 



(5.28) 



d 2 



'ab-^z+Vu,-^:. The corresponding boundary conditions read /x a (|z| 



where V ab = D 

oo) = 0. The matrix valued Greens function G a b is defined by ^2^ =1 (—k ac V ac + 
T ac )G c b(z) = 5 a b5(z). In the transport equations ( |5.23| ) position dependent rates 
are present, e.g. T m . They typically vanish in the symmetric phase and be- 
come maximal in the broken phase. In order to keep the problem analytically 
tractable we simply model the position dependence of these rates by step func- 
tions, i.e. T ah (z) = T +ab Q(z) + T_ ab Q(-z). 



6 Numerical results 

In this section we present our numerical results for the baryon asymmetry in the 
MSSM, which has to be compared with the observational value 2 — 7 x 10~ n [p4}| . 
We model the Higgs field H\(z) in the bubble wall by a kink with width L w and 
take H 2 (z) = ta,nf3(z)H 1 (z), where tan/3 varies in the wall. In all our evaluations 
we assume v c = 120 GeV and tan/3 (T c ) = 3. For the critical temperature we take 
T c = 90 GeV. 

In Section 5.3 we emphasized the impact of the squark spectrum on the baryon 
asymmetry. A strong phase transition requires the right-handed stop to be light. 
In the transport equations we simply treat it as a massless particle, i.e. kx = 18. 
We assume all other sfermions to be heavy compared to T c , i.e. fcg 3 = ks = 
9. From Eq. ( p.24| ) we find for the effective statistical factors in the transport 
equation A = 126/13 and B = 23/13, and from Eq. ( gTggj ) C = 5/13. This 



non-universal squark spectrum has the additional virtue of relaxing the strong 
sphaleron suppression of the baryon asymmetry [ 55 1 . The Higgs fields are treated 
as massless as well, i.e kn = 6. 
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Figure 9: Hi ± H2 contributions to the baryon asymmetry dependent on v w for 
different values of the wall thickness L w = 20/T, 15/T, 10/T (from below) and 
= |M 2 | = 150 GeV, arg(/iM 2 ) = n/2 and 5/5 = 0.01. r\ is given in units of 
2 x 10~ n (observational bound). 



For the diffusion constants of quarks and Higgses we take Dh = 110/T and 
D q = 6/T. The rates T v a and T y ^ which involve heavy squarks are Boltzmann 
suppressed and we set them to zero in the following p3| . Furthermore we take 
T y = 0.015T, T hf = 0.016T, T m = 0.05T9(-z), T Hl = T H2 = 0MT6(-z). The 
weak sphaleron rate is r„, s = 20a^,T. 

In the following we separately present the Hi—H 2 and Hi + H 2 contributions to 
the baryon asymmetry, which originate from the flavor and helicity contributions 
to the chargino dispersion relation, respectively. We use the dispersion relations 
in the canonic momentum. In the formulation with the kinetic momentum the 



Hi — H 2 part vanishes [18]. However, a Hi — H 2 contribution has been found 
in real time Green's function treatments of the chargino current, most recently 



in ref. |[23|| . It would we very interesting to understand more clearly the relation 
between our Hi — H 2 source and the one found in that approach. 

In fig. |9] we summarize our results for the baryon asymmetry generated during 
the phase transition. To maximize the result we take \fi\ = \M 2 \ and maximal 
CP violation aig(M 2 fi) = tt/2. For smaller phases the result simply scales with 
sin(arg(M 2/ u)). Moreover, the Hi — H 2 result is proportional to the change in 
the Higgs vev ratio in the bubble wall, which we take to be 5/3 = 0.01. We work 
with fx = 150 and tan/3 = 3. The mass of the lightest chargino eigenstate is then 
about 110 GeV, in agreement with the experimental constraints. 

In our evaluations we vary the wall velocity in the interval 10~ 5 < v w < 0.5. 
For small values of v w we deal with a quasi equilibrium situation, and the baryon 
asymmetry goes to zero. Also large wall velocities suppress 77, since transport 
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in front of the bubble wall becomes inefficient. From Fig. |9] we observe that 
baryogenesis is most efficient if v w is a few times 10~ 2 . This behavior was already 



found in refs. [18,23]. Most interesting, calculations of the wall velocity in the 



MSSM lead to v w indeed in that range [ 1^, p3 1 . The detailed calculations of 
ref. |25] lead toQ v w = 2 — 3 • 10~ 2 , very near to the maximum for 77. 

We also present results for different wall widths L w = 10/T, 15/T, 20/T. As the 
derivatives with respect to the z coordinate in the source term suggest, thinner 
walls induce a larger baryon asymmetry. However, this effect is much more 
pronounced for the H\ — H 2 contribution. 

For the Hi + H 2 contribution we find a baryon asymmetry which in the most 
favorable case is about 7 times larger than the observed value. Using kinetic 
momentum in the averaging, one obtains maximally (in the massless case) a 
further factor 2. As result, the assumption of maximal CP violation can only 
mildly be relaxed. Large complex phases are necessary, which are only compatible 
with EDM experiments if accidental cancellations happen or the first and second 
generation squarks have TeV scale masses. The H\ — H2 contribution to rj is 
an order of magnitude larger, and smaller phases of the order 10~ 2 would be 
sufficient. 



Let us finally compare our results with those of refs. [|11|[23|], which also started 
from the transport equations ( |5.23|) . The authors of ref. [18] also use the semi- 
classical method but consider only the if 1 + H2 contribution, and take somewhat 
different values for the interaction rates and diffusion constants. Our results 
are qualitatively similar .0 The dependence of the baryon asymmetry on v w shows 
small deviations. This may be due the fact, that in contrast to ref. [18] we exactly 
calculate the Greens function associated with the transport equation (|5.23|). 



Ref. [23] is based on quantum transport equations. In this derivation of the 
chargino source terms also a Hi — H 2 contribution appears. The results in am- 
plitude and velocity dependence are very similar to those presented here for the 
H\ — H 2 source. The small deviations in the regime of small v w may again be 
due to an approximation in the computation of the Greens function if ref. p3| . It 



remains to be seen if the Hi — H 2 contributions in the approach of refs. [21-23 



and in the semi-classical method are related, and if they are indeed artificial as 
suggested by the use of kinetic momentum. 



4 Potential changes in v w due to reheating [g6| were neglected. 

5 taking into account a numerical normalization error corrected in a forthcoming revised 
version of (lal. 
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7 Conclusions 



In the first part we have given in some detail the procedure of solving numerically 
for the wall profile in the case of more then one Higgs field for a given effective 



potential. This extends a previous description |13|]. Within this method we can 
also systematically search for CP violating solutions. It turns out - and this is 
confirmed in ref. [0] - that in the MSSM there is no indication of a spontaneous 
CP violation. However, in the larger parameter space of the NMSSM |2S| such 
solutions can be found in 13]. In the MSSM we remain with the possibility of 
explicit CP violating interactions, in particular of the charginos. 

In the second part we treat transport around the stationary proceeding wall 
in the quasi-classical approximation. The solution of the set of Dirac equations 
for the charginos in the WKB approximation shows a split in the dispersion re- 
lations between particles and anti-particles at order h because of CP violation. 
We have written the full dispersion relations according to ref. [[nj for canoni- 
cal momenta, but also presented the relations for kinetic momenta. In ref, 



arguments have been presented that one should use the kinetic momentum in 
Boltzmann transport equations. This would destroy all source terms of CP vio- 
lation antisymmetric in the two Higgses which have been considered to be most 
important in older work on the MSSM. Also in order to compare with work in 
the spirit of quantum Boltzmann equations [pi"1-|2"3"|j we have kept open for both 
versions. 

Of course one also has to test the two different types of momentum in averaging 
procedures, but this just gives a factor (about 2). Indeed, because tan/3 only varies 
very weakly in the wall, the contribution of the symmetric Higgs combination can 
not be neglected anyway. However, if the antisymmetric source is absent, rather 
large explicit CP violation is needed to obtain the observed baryon asymmetry 
even for favorable parameters of the MSSM. Thus the discussion of consistency 
with experimental EDM bounds due to accidental cancellations becomes relevant. 
On the other side this might point to an extension of the MSSM, e.g. to some 
version with a singlet. It is encouraging that we at least obtained a baryon 
asymmetry which is of the right order of magnitude compared to observations. 
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